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First-principles calculations are performed to investigate f magnetism in A 2 Ti 2 07 (A=Eu, Gd, 

Tb, Dy, Ho, Er, Yb) magnetic pyrochlore oxides. The Hubbard U parameter and the relativistic spin 
orbit correction is applied for more accurate description of the electronic structure of the systems. 

It is argued that the main obstacle for first-principles study of these systems is the multi-minima 
solutions of their electronic configuration. Among the studied pyrochlores, Gd 2 Ti 2 07 shows the least 
multi-minima problem. The crystal electric held theory is applied for phenomenological comparison 
of the calculated spin and orbital moments with the experimental data. 


I. INTRODUCTION 

Magnetic pyrochlore oxides [l] with chemical formula 
A2B2O7 have rich physics and exotic magnetic properties 
(such as magnetricityQ) caused by geometrical frustra¬ 
tion. In these materials A and B are usually trivalent 
rare-earth and travalent transition metal ions, respec¬ 
tively, which form a network of corner sharing tetrahedral 
on the fee Bravais lattice (Fig. []}. This geometrical fea¬ 
ture is known to be the origin of the magnetic frustration 
of the system. Among these materials, A 2 Ti 2 07 (A=Eu, 
Gd, Tb, Dy, Ho, Er, Yb) are the most studied compounds 
[3], and exhibit various magnetic phenomena including 
spin ice behavior in Ho2Ti207 and Dy 2 Ti 2 07 jjj, spin 
liquid behaviour in Tb2Ti2C>7 [U, [p and Yb2Ti2C>7 3, 
and order by disorder phenomena in Er2Ti207 [l|, LZI, ■ 
Gd 2 Ti 2 07, as a controversial case, is assumed to be a 
classical Heisenberg antiferromagnet [oj. 

Theoretical investigations on magnetic pyrochlores are 
mainly based on model Hamiltonians, involving Heisen¬ 
berg exchange, dipole-dipole interactions, and single ion 
anisotropy, which are parametrized by using experimen- 



FIG. 1: The structure of pyr-A 2 Ti 2 07 . Left: positions of 
magnetic A 3+ ions (red color) and Right: positions of Ti 4+ 
ions (blue color) on the the cubic pyrochlore lattice. The O -2 
ions (not shown) are located at the bent (in gray) where bonds 
merge. 
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tal data pj. The lack of modern ab initio calculations 
on these materials prevents microscopic understanding of 
these magnetic model Hamiltonians. To our knowledge, 
ab initio calculations on magnetic pyrochlore oxides are 
limited to nonmagnetic properti es |ToL 0 or those com¬ 
pounds with no active f electrons [l2| • It is due to the fact 
that the partially filled 4f orbital of the rare earth ions 
involves strongly correlated and localized electrons which 
cause serious difficulties in finding the true ground state 
of the system within density functional theory (DFT) 
computations [13]. It is well understood that conven¬ 
tional local functionals give rise to wrong ground states 
for these strongly correlated systems. The usual solu¬ 
tion for correct description of the coulomb interaction 
between 4f electrons is applying orbital dependent ap¬ 
proaches including hybrid functionals [I3| and the Hub¬ 
bard based DFT+U technique. An important challenge 
in first-principles calculation of the systems with 4f elec¬ 
trons within orbital dependent functionals is sensitivity 
of the final results to the initial electronic configuration of 
the 4f shell. Moreover, close energy local minima are se¬ 
rious obstacles to achieve convergency in first-principles 
calculation of the 4f electron systems. 

Our specific aim in this work is employing DFT+U ap¬ 
proach for first-principles calculation and investigation of 
f magnetism in A 2 Ti 2 07 (A=Eu, Gd, Tb, Dy, Ho, Er, 
Yb) compounds. After brief explanation of our com¬ 
putational method, the obtained results are presented in 
section III. In section IV, our conclusions are presented. 

II. COMPUTATIONAL METHOD 

Our electronic structure calculations are performed by 
using Fleur [l5| computer code, which solves Kohn-Sham 
single particle equations by using the full-potential lin¬ 
earized augmented plane wave (FP-LAPW) method IHt . 
The muffin-tin radii of A, Ti, and O atoms were set 
to 2.75, 2.0, and 1.35 bohr, respectively. The Perdew- 
Burke-Ernzerhof (PBE) formulation of the generalized 
gradient functional (GGA) is used in this work E3- We 
used the GGA+U approximation for better description 
of the Coulomb interaction between f electrons. The 
adapted values of U for A 2 Ti 2 07 are given in table HI 
while the value of on-site Hund’s exchange, J h , is set to 
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TABLE I: Hubbard U parameters used in this work (taken 
from Ref. 0), calculated band gap (this work), lattice con¬ 
stant (a), and structural internal parameter of oxygen ( x) 
0 in the studied titanate pyrochlore oxides. The available 
experimental values of the band gaps are given in the paren¬ 
thesis. 



U (eV) 

gap (eV) 

a (A) 

X 

Eu. 2 Ti. 2 O 7 

10.14 

2.1 (2.5 [20]) 

10.194 

0.327 

Gd2Ti207 

11.48 

3.0 (3.2 [211) 

10.186 

0.326 

Tb2Ti207 

5.00 

2.4 

10.159 

0.328 

Dy2Ti207 

5.68 

1.9 (2.4 [22]) 

10.124 

0.328 

Ho2Ti207 

6.82 

- (3.2 [23]) 

10.104 

0.329 

Er2Ti207 

6.80 

3.0 

10.071 

0.328 

Yb 2 Ti 2 0 7 

6.00 

0.6 

10.325 

0.331 


1 eV. The spin-orbit coupling (SOC) is very important in 
the systems with heavy 4f elements, hence we considered 
this relativistic interaction in our calculations. 

A 4 x 4 x 4 Monkhorst-Pack k mesh was used for Bril- 
louin zone integration. The lanthanide 4/ electrons were 
treated as valance states while their 5s 2 5p 6 electrons as 
well as the Ti 3s 2 3p 6 electrons were considered as semi¬ 
core states. The cut-off of wave function expansion in the 
interstitial region was set to 3.8 (a.u) -1 . We adapted ex¬ 
perimental lattice parameters of A 2 TUO 7 , given in Table 
HI for our calculations. The space group of A 2 Ti 2 07 is 
Fd3m, with the following Wyckoff positions for A and Ti 
atoms, receptively; 16d (1/2,1/2,1/2), Ti: 16c (0,0,0). 
There are two Wyckoff positions for oxygen; O at 48f 
(x, 1/8,1/8), and O' at 8 b (3/8, 3/8, 3/8). The primi¬ 
tive cell of A 2 U 2 O 7 , contains 22 atoms and the internal 
parameter x was taken from experiment (Table HI) . 

Unfortunately the calculations for E^T^CUand 
Ho 2 Ti 2 07 were unstable and we could not achieve their 
converged electronic structure within GGA+U+SOC. 
Hence in the case of Er 2 Ti 2 07 , we report our previous 
converged results [24j within LDA+U+SOC and muffin- 
tin radii 2.5, 2, 1.5 bohr for Er, Ti, and O, respectively. 
While for Holmium Titanate, we could converge its elec¬ 
tronic structure within GGA+U to determine its band 
gap and spin moment. 

III. RESULTS 

As it was mentioned, the multi minima problem is 
a serious challenge in first-principles investigation of 
4f electron systems. In order to see the existence of 
this problem, we performed some test calculations on 
Tb 2 Ti 2 07 within GGA+U, by starting from three dif¬ 
ferent initial states. In these preliminary calculations, 
we omitted the SOC interaction, while other parts of our 
study involves this relativistic correction. At the end of 
calculations, we reached to three different electronic band 
structures, presented in Fig. [2] The top band structure 
is obtained by using the full symmetry (Fd3m) of the 
system during the computations, while the middle band 
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FIG. 2: Obtained three different band structures of 

Tb 2 Ti 2 07 within GGA+U. Please see text for more details. 
The Fermi energies are set to zero. 


structure is calculated without using any symmetry. It 
is clearly visible that in the resulted band structure af¬ 
ter low symmetry calculations, the narrow occupied and 
unoccupied / bands are slightly separated from the va¬ 
lence and conduction bands. The bottom band structure 
was calculated by using the full symmetry of the struc¬ 
ture and the converged metallic electronic structure of 
the system within GGA as the starting point of GGA+U 
calculation. In this way, the GGA+U calculation is con¬ 
verged to a metallic electronic structure. 

As it was mentioned, GGA predicts a wrong metal¬ 
lic state for most of the 4f compounds, while GGA+U 
opens a gap between f states and make A 2 Ti 207 insula¬ 
tors. The calculated band gap of these compounds within 
GGA+U+SOC are given in TableQ] along with the avail¬ 
able experimental data. The good agreement observed 
between computed and measured band gaps confirms the 
selected U parameters for these materials. The obtained 
orbital resolved Density of State (DOS) of the studied 
systems are shown in Fig. [3j As it was mentioned in 
Computational method, there are two kinds of oxygen 
atoms in the systems. The first kind of O atoms con¬ 
nect the neighboring Ti 4+ ions displayed in Fig. |T] and 
hence, a good hybridization happens between d-orbital 
of Ti and p-orbital of these O atoms (Fig. [3]). The sec- 
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FIG. 3: Calculated partial density of states (DOS) of the 
studied titanate pyrochlore A 2 Ti 207 (A=Eu, Gd, Tb, Dy, 
Yb) by GGA+U+SOC. The Fermi energies are set to zero. 


ond kind of O atoms (O') are inside tetrahedrons of A 3+ 
ions, hence their p orbital DOS are compared with the 
f orbital DOS of A atoms in Fig. [3] The hybridization 
between these orbitals is clearly week, which is likely due 
to the localized nature of f states. However, in the case of 
Tb, Dy, and Yb based pyrochlores, the contribution of f 
states in the valence shell is high, showing the importance 
of valence treatment of f electrons for ab initio calcula¬ 
tion of these systems. It is seen that among A 2 Ti 207 , 
Gd 2 Ti 2 07 exhibits the lowest contribution of f states in 
the valence shell, which is due to the half-filling of f or¬ 
bital in Gd. As a result of that, there is only one set of 
possible occupation numbers in Gd 2 Ti 207 and this helps 
the ab initio calculations to easily converge to its global 
minima. 

The calculated spin, orbital, and total magnetic mo¬ 


ments of the studied A 2 Ti 207 pyrochlores obtained from 
collinear GGA+U+SOC calculations are given in Ta¬ 
ble m These moments are calculated inside Muffin-tin 
spheres. The spin quantization axis is set to the z direc¬ 
tion. The experimental values of total moments are also 
reported for comparison. The calculated effective mag¬ 
netic moment of Eu 2 Ti 207 is about 4.9 hb, which shows 
large deviation from the experimental value of 0.6/is [25j . 
In contrast, the calculated energy gap of this system is 
close to the experimental value (Table Q]). In the case 
of Gd 2 Ti 207 we got a magnetic moment of 7.0 ^b, in 
good agreement with the measured value of 7.7hb IS!])- 
The Gd 3+ ion with a valence shell of (4 f 7 ) has a spin of 
S = 7/2 with no orbital magnetic moment, giving rise 
to an effective moment of g x S(S + 1) = 7.9 (g=2), 
in good agreement with the experimental moment. The 
GGA band gap of Gd 2 Ti 2 07 is found to be about 2.7 
eV, in agreement with Xiao et al. [26|. This compound 
is the only titanate pyrochlore with a band gap within 
GGA. This observation provides further evidence for the 
fact that this system easily converges to its true ground 
state during ab initio calculations. It should be admitted 
that the true ground state of this system, as well as other 
magnetic pyrochlores, may be a non-collinear antiferro¬ 
magnetic state. Our preliminary calculations confirm the 
non-collinear ground state of Gd 2 Ti 2 07 , however accu¬ 
rate determination of this ground state requires compre¬ 
hensive ab initio calculations which is out of the scope of 
this paper. It should be noted that non-collinear calcula¬ 
tions for other pyrochlores are much more complicated. 

In Tb 2 Ti 2 07 , GGA+U+SOC gives a total moment of 
about 7.2/rs, which is larger than the measured value of 
5.1 hb 0- On the other hand, the calculated total mo¬ 
ment of Dy 2 Ti 2 07 (8.2 ^b) is lower than the experimental 
moment of 10 ^b ;3J- la the case of Yb 2 Ti 2 07 , our ab 
initio calculation predicts a moment of 2.4 hb which is 
about twice the experimental value of about 1.1 /is mea¬ 
sured by single-crystal neutron diffraction method be¬ 
low 0.2 K [6]. As it was mentioned in Computational 
method, we could not converge our GGA+U+SOC calcu¬ 
lations for Er 2 Ti 2 07 and Ho 2 Ti 2 07 . The total moment 
of Erbium Titanate within LDA+U+SOC was found to 
be 8.8 [Ab, which is about 5.7 fiB bigger than the experi¬ 
mental value [§]. The spin moment of Ho 2 Ti 2 07 within 
GGA+U was found to be about 3.8 /x_b 

For more accurate comparison of the experimental and 
computational values of magnetic moments, we used the 
crystal electric field theory to decompose experimental 
moments into spin and orbital contributions. Apply¬ 
ing crystal electric field theory to the experimental neu¬ 
tron spectra, the magnetic ground-state wave function of 
A 2 Ti 2 07 (A=Eu, Gd, Tb, Dy, Ho, Er, Yb) pyrochlore 
compounds, except Gd 2 Ti 2 07 and Eu 2 Ti 2 07 are deter¬ 


mined in terms of | j,rrij) basis set 


where j and 


mj are total angular momentum and magnetic quantum 
number, respectively. We used Clebsch-Gordan coeffi¬ 
cients to decompose \j,nij) to \mi,m s ) states, where mi 
and m s are orbital and spin magnetic quantum num- 
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TABLE II: Magnetic properties of titanate pyrochlore, L,S,J: orbital, spin, and total angular momentum, ns, hl (mb) : 
Calculated spin and orbital magnetic moments, The values in the parenthesis are estimated from crystal field analysis, fitot 

M : Calculated total magnetic moment, The values in the parenthesis are available measured data, collected from references 

MMMM 



L 

S 

J 

Ms 

I^L 

fJ'tot 

Eu. 2 Ti. 2 O 7 

3 

3 

0 

6.1 (-) 

- 1.1 (-) 

4.9 (0.6) 

Gd2Ti207 

0 

7/2 

7/2 

7.0 (-) 

0.0 (-) 

7.0 (7.7) 

Tb2Ti207 

3 

3 

6 

6.0 (3.3) 

1.2 (1.7) 

7.2 (5.1) 

Dy2Ti207 

5 

5/2 

15/2 

5.0 (4.9) 

3.3 (4.9) 

8.2 (10.0) 

Ho2Ti207 

6 

2 

8 

3.8 (3.9) 

- (5-8) 

- (10.0) 

Er2Ti207 

6 

3/2 

15/2 

3.0 (1.2) 

5.9 (1.9) 

8.8 (3.2) 

Yb2Ti207 

3 

1/2 

7/2 

0.9 (0.3) 

1.5 (0.8) 

2.4 (1.1) 


bers, respectively. The resulting magnetic ground-state 
wave functions in terms of \mi,m s ) states are given in 
appendix [A] and the obtained orbital and spin moments 
are presented in Table |TT] The observed consistency be¬ 
tween crystal field spin and orbital moments and mea¬ 
sured total moments confirms the obtained magnetic 
ground state wave functions. It is seen that the calcu¬ 
lated orbital moment of Tb 2 Ti 2 C >7 and Dy 2 Ti 2 C >7 are 
about 30% smaller than the crystal field orbital moment, 
while the calculated orbital moment of Er 2 Ti 207 and 
Yb 2 Ti 2 C >7 are significantly higher than the crystal field 
values. On the other hand, the calculated spin mo¬ 
ments show generally more agreement with the corre¬ 
sponding crystal field values. The obtained spin moments 
of Dy 2 Ti 2 07 and Ho 2 Ti 2 07 are very close to the crystal 
field data while in the case of Tb 2 Ti 207 , Ei’ 2 Ti 207 , and 
Yb 2 Ti 207 significant difference is visible between first- 
principles and crystal field spin moments. 


IV. CONCLUSIONS 

The electronic structure and magnetic properties of 
A 2 Ti 2 07 (A=Eu, Gd, Tb, Dy, Ho, Er, Yb) magnetic 
pyrochlores were investigated by using full-potential den¬ 
sity functional calculations within GGA+U scheme and 
relativistic spin orbit coupling (SOC). The calculated 
band gaps are in good agreement with the available mea¬ 
sured values. We used Tb 2 Ti 2 07 compound to show the 
possible multi-minima solutions of the electronic struc¬ 
ture of these magnetic pyrochlores. It was argued that 
the half filling of the 4f shell of Gd significantly de¬ 
creases the contribution of f electrons in the valence shell 
of Gd 2 Ti 2 07 and hence facilitates the electronic struc¬ 
ture calculation of this system to converge to its global 
minima. The calculated spin, orbital, and total mag¬ 
netic moments of the systems within GGA+U+SOC cal¬ 
culations were presented and compared with the avail¬ 
able measured total magnetic moments. The best agree¬ 
ment was observed in Gd 2 Ti 207 , while Eu 2 Ti 207 and 
Ei^T^O 7 compounds showed the highest deviation from 
experiment. The phenomenological magnetic ground 
state wave function of Tb 2 Ti 2 07 , Dy 2 Ti 2 07 , Ho 2 Ti 2 07 , 


Er 2 Ti 2©7 , and Yb 2 Ti 207 compounds were determined 
in the framework of crystal electric field theory for de¬ 
composition of the measured total magnetic moments 
into spin and orbital contributions. Comparing the cal¬ 
culated and the crystal field spin and orbital moments 
shows more accuracy of the calculated spin moments. 
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Appendix A: Crystal field analysis 

The magnetic ground-state wave function of A 2 TUO 7 
(A=Eu, Gd, Tb, Dy, Yb) pyrochlore compounds in terms 
of \j,rrij) states are as follows [27 ]: 

\(j> T g b ) = +0.96|±4)±0.13|±l) 

—0.12| =F 2) =F 0.23| =F 5> 

1 ^ Q 

K y ) = +0.981I ±y)± 0.1901 ±-> 

+0.022|±^)± 0.037| 

+0.005| + ^)±0.001| + y) 

| <t>f°) = -0.979| ± 8) ± 0.189| ± 5) 

-0.014| ±2) ±0.070| + 1) 

-0.031| + 4) ± 0.005| + 7) 

IO = + 0.376| ± 7 -) + 0.9221 ± ±) 

- 0.093| + |> 
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\<t>g r ) = + 0-471 ±y}± 0.42| ±1) 

- 0.57| ±1) T 0.24|t^) 

+ 0.47|=F y) 

We used Clebsch-Gordan coefficients to represent 
these magnetic ground-state wave functions in terms of 
| mi,m s ) states: 

I <% b ) = + 0.46| - 3, -1) + 0.711 - 2, -2) 

+ 0.46| - 1,-3) - 0.04| + 1,-2) 

— 0.08| 0,-1)-0.08|-1, 0) 

— 0.04| — 2, +1) — 0.011 + 2, —3) 
-0.01|-3, 2) — 0.021 3,-1) 

— 0.06| 2, 0) -0.08| 1, 1) 

— 0.06| 0, 2) — 0.021 - 1, 3) 

+ 0.16| 3, 2) + 0.161 2, 3) 

| (j)f y ) = + 0.9800| ± 5, ±|) ± 0.0970| ± 2, ±|) 

± 0.1330| ± 3, ±|) ± 0.0890| ± 4, ±i) 

± 0.0280| ± 5, t\) T 0.00101 T 4, ±|) 

± 0.0080| T 3, ±|) ± 0.0180| T 2, ±^) 

± 0.0240| T 1,7^) ±0.0180| 0,t|) 

± 0.0070| ±1,t|)- 0.0045| Tl,±|) 

-0.01101 0, ±|> - 0.01421 ± 1, ±i) 

- 0.0107| ± 2, t\) - 0.0046| ± 3, t|) 

- 0.0098| ± 4, t|) + 0.0007| T 5, ±^) 

+ 0.0023| T 4, t\) ± 0.0035| T 3, t|) 

+ 0.0025| T 2, t|)± 0-00101 T 5, t\) 

\<f>£°) = - 0.9790| ± 6, ±2) ± 0.1184| ± 3, ±2) 

± 0.1290| ± 4, ±1) ± 0.0670| ± 5, 0) 

± 0.0160| ± 6, =Fl) - 0.0047| 0, ±2) 

- 0.0088| ± 1, ±1) - 0.0080| ± 2, 0) 

- 0.0046| ± 3, =Fl) - 0.0012| ± 4, =p2) 

± 0.0098| T 3, ±2) ± 0.0300| T 2, ±1) 

± 0.0450| =F 1, 0) ± 0.0400| 0, =f1) 

± 0.0180| ± 1, t2) - 0.0007| T 6, ±2) 

- 0.0050| T 5, ±1) - 0.0210| T 3, ±1) 

- 0.0160| T 2, ±2) - 0.0140| T 4, 0) 

± 0.0025| T 6, ±1) ± 0.0043| T 5, ±2) 


lO = + 0.37| ± 3, ±i)± 0.70| 0,±i) 

± 0.60| ±1,T^)-0.03| T 3, ±i) 
— 0.091 =F 2, =pi) 


= + 0-21 1 ± 6, ±\) ± 0.42| ± 5, ±|) 
±0.04|±5,t^)±0.16|±4,t^) 
± 0.31| ± 3, ±|) ± 0.25| ± 2, ±|) 

-0.16|±2,t|)-0.35|±1,t^) 

- 0.38| 0, ±|) - 0.20| =F 1, |> 

T0.12|Tl,T|)T0.17|T2,T^) 

T 0.111 T 3, ±i) =F 0.04| =F 4, ±^) 
+ 0.37| T 4, t|) + 0.28| =F 5, 
±0.08|t6,±^) 


For calculating the spin and orbital magnetic moments 
from the above wave functions, one should take into ac¬ 
count that the Tb 2 Ti 2 07 , Dy 2 Ti 2 07 , Ho 2 Ti 2 07 , and 
Yb 2 Ti 2 C> 7 are Ising like systems while Er 2 Ti 2 C >7 adopts 
XY model magnetism. For Ising like systems, gi(L z ) and 
g s (S z ) are calculated as orbital and spin moments of the 
magnetic atom of the systems, while in XY-like systems, 
(S±) » {S z ) and (L±) > (L z ), therefore gi{L±) (L± = 
L x ± iL y ,gi = 1) and g s (S±) (S± = L x ± iL y ,g s = 2) are 
calculated as orbital and spin moments, respectively. 
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